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ABSTRACT 

Method A method of generating a hybrid grid of a heterogeneous formation 
crossed by one or more geometric discontinuities such as, for example, an underground 
formation where one or more wells have been drilled, or a fractured formation, by 
combining structured grids and non-structured grids in order to carry out simulations in 
accordance with a defined numerical pattern is disclosed. 

__ Hybrid gridding in an application to a medium crooood by wclla for example i s 

c^cntially performed by associating a first structured grid (Gl) for gridding of the 
heterogeneous medium rnr . pnr.ting the considering discontinuities thereof with second 
structured, radial type grids (G2) for gridding of a zone around each pipe or well, which 
allows te-better pirtimlnr consideration of constraints linked with flows in this 

thezone. m order to connect the first grid of the medium and the second well grids, non- 
structured transition grids (G3) are interposed there between thorn and th e ^Apower 
diagram technique is used-therefefe, which is particularly advantageous in that it allows 

nrirnrrinHy ""^nr t appro priate connection of n on-regular structured grids._An 

application is 



Applications hydrocarbon reservoir simulators for oxample simulati 
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BACKGROUND OF THE INVENTION 
nT 7| p r>T7 t i te rMWFMTTn W Field of the Invention 
The present invention relates to a method of generating a hybrid grid of a 
heterogeneous formation crossed by one or more geometric discontinuities in order for 
example to carry out simulations for example . 

The method is more particularly applied to formation of a grid suited to an 
underground reservoir crossed by one or more wells, or by fractures or faults, in order to 
model displacements of fluids such as hydrocarbons. 

BACKGROUND ™ ™p rMVFNTTONDescription of the Prior Art 
Grid generation is a crucial element for the new generation of reservoir simulators. 
Grids allow to describe the geometry of the geologic structure studied by means of a 
representation in discrete elements wherein simulation is performed according to a suitable 
numerical pattern. Better comprehension of physical phenomena requires 3D simulation of 
the multiphase flows in increasingly complex geologic structures, in the vicinity of several 
types of singularities such as stratifications, faults, pinchouts, channels and complex wells. 
All this complexity has to be taken into account first by the grid which has to reproduce as 
accurately as possible the geologic information in its heterogeneous nature. 

Grid modelling has made great advances during the past few years in other fields such 
as aeronautics, combustion in engines, structure mechanics, etc. However, the gridding 
techniques used in the other fields cannot be applied as-they-arfr-to th^petroleum sphere 
ap plications b ecause the professional constraints are not the same. For example, in reservoir 
simulation, the numerical patterns are constructed from control volumes in order to better 
respect the mass conservation in the case of transport equations of hyperbolic nature. The 



grid must be a Mnnir ™nt„r«H ^ "block-centered" type grid, i^thatjs the nodes must be 
inside each layer and the boundaries of each block must follow the interface between the 
layers. Now, if this constraint was not taken into account, the nodes would naturally be 
placed along the faults and along the stratification boundaries. The consequence of this would 

5 be that these interfaces would pass through the control volume used. The saturation, constant 
in the control volume, could not rcopoct then the diocontinuity consider discontinuities a nd 
the results would not be accurate. It is therefore necessary to develop new techniques that are 
better suited to the-petroleum sphere-a pplication r equirements. 

Cartesian grids, which are commonly used in current commercial simulators, are 

10 unsuited for solving these new problems posed by the development of petroleum reservoirs. 
Cartesian grids, based on parallelepipedic elements, do not allow representation of such 
complex geometries. 

There is a well-known method of generating structured 3D hexahedral grids of CPG 
(Corner-Point-Geometry) type which respects the geometry of the bodies. It is described in 
1 5 French p atent FR-2,747,490 f US -patent 5,844,564) filed by the appMeafl^assignee_and also 

in the following publication: 

Bennis Ch. Et al. ^ One More Stop in Gocad Stratigraphic Grid Goncration^ _^One 
More Step in Gocad Stratigraphic Grid Generation : Taking into Account Faults and 
Pinchouts^ ; SPE 35526, Stavanger, 1996. 
20 This grid type is more flexible than cartesian grids because it consists of any 

hexahedral elements that can be degenerated. It strictly respeets xonsiders t he horizons, the 
faults and it allows tn mpmrrmt re presentation of certain unconformities such as pinchouts 
because its construction is based on these elements. However, this type of grid does not allow 
te^elve-solutionof all the-geometric complexities such as, for example, circular radial grids 



around complex wells. It is possible to form separately the grid of the reservoir and the grids 
around the wells but it is difficult to represent several objects in the same CPG type reservoir 
grid because of connection problems linked with the structured nature of the grid. 

Another approach is also known where 3D grids arepnly based on tetrahedral 
5 Delaunay elements, with a circular radial refinement around the wells, are-bejng 

automatically generated. The advantage of such an approach is that it is entirely automatic 
and does not practically net-require the user's attention. However, this method has drawbacks 
which make the results obtained difficult to use: 

- there are on average five times as many grid cells as in a CPG type grid for the same 
10 structure, which is very disadvantageous for simulation calculations, 

- unlike the structured grids which are easy to-^4sualiz©display, to explore from the inside and 
to locally modify interactively, it is very difficult and sometimes impossible to properly 
control the tetrahedral grids because of their size and especially because of their non- 
structured nature. This poses problems for validating the grid from a geometric point of view 

1 5 as well as for understanding and validating the result of a simulation on this type of grid. 

Other approaches are also well-known, which allow to g e n e rat e generation of g rids, 
notably grids based on control volumes generated from a triangulation associated with 
techniques of aggregation of the triangles (or tetrahedrons) into quadrangles allowing the 
number of grid cells to be reduced. Although promising results were obtained with these new 
20 grids, precise representation of the geologic complexity of reservoirs and wells remains a 

subject for research and development. Despite their hybrid-aspeetiiature, they remain entirely 
unstructured and would therefore be very difficult to manage and to handle in real 3D. 
Furthermore, taking account of real 3D faults and deviated wells would greatly increase this 
difficulty. 



Frenc^Patent-EatenLapplication FR-99/1 5,120, filed by the appHeant-assignee, 
describes a method of generating a 3D hybrid grid on a heterogeneous medium, comprising 
using a first structured grid, of CPG type for example, and structured radial grids around well 
or pipe trajectories drawn or imported m-into the model in order to better take into account 
5 the particular constraints linked with the flows in the vicinity of these wells. They- These grids 
are combined by including the radial grids around the-ei-around each well in the global 
reservoir grid after forming therein cavities that are large enough to allow formation of «eft- 
Ktriictured unstructured transition grids. 

There are various well-known techniques for forming non otructuredunstructured 
1 0 grids. Tb«v -These techniques can be based on canonical polyhedrons (tetrahedrons, 

pentahedrons, pyramids, etc.) according to the numerical pattern beingused, and entirely 3D 
solutions applicable to these grid types are known. 

The method according to the invention is intended for 2.5D generation of transition 
grids based on Voronoi type or similar polyhedrons in order to apply control volume type 
1 5 numerical patterns, which are reduced to a 2D problem provided that all the layers of the well 
grids and of the reservoir grid can be projected vertically in a horizontal plane so as to form 
identical grid cells. 

Such a grid must meet certain conditions: 

the segment connecting the eenftes -centers o r sites of two noighbouring neighboring 
20 cells must be orthogonal to the side common to the two cells, 

two nnighhnurinp neighboring cells must entirely share the side that connects them, 
fc fethat is a cell can 5 be connected to only one other cell by the same side, and 

the cells must be convex. 
In 2D, these constraints impose that: 



the edges of the boundary polygons must not be modified (each one must correspond 
to an edge of a cell created, without subdivision of the edge), and the cells which are c reated 
must be convex, 

each cell must have a eentre- center such that the straight line connecting the contreo 
5 centers of two adjacent cells is perpendicular to the edge shared by the two cells, and 
the centres- centers of the cells must be inside their cell. 
Other entirely non-structured approaches are also well-known, which allow to g e n e rate 
generation of g rids, notably grids based on control volumes generated from a triangulation, 
associated with techniques of aggregation of the triangles (or tetrahedrons) into regular 
1 0 quadrangles allowing the number of grid cells to be reduced. 

Voronoi diagrams are formed, as it is well-known, by determining ar6und a set of 
sites P distributed in a space regions consisting of series of points which are closer to each 
site than to any other site P of the set. This construction technique is described in detail 
notably in the following documents: 
15 - Aurenhammer F., 1991 : Voronoi diagrams "Diagrams : arA_Survey of Fundamental 
Geometric Data Structure,!! ACM Comput. Sun^, 23, 345-405, 

- Fortune S., 1992 : Voronoi Diagrams and Delaunay Triangulations, pages 225-265 of D.Z. 
Du & F.K. Hwang (eds), Computing in Euclidean Geometry, 2nd edn., Lecture 
Notes Series on Computing, vol.4, Singapore, World Scientific. 
20 The cells of a Voronoi diagram meet the following constraints : the orthogonality 

between the 6entfe&^enters_and the adjacent edge is respected and the cells are convex. It is 
therefore natural to propose filling the cavity with Voronoi cells while trying to respeet 
consider the edge geometric constraints. However, if a Voronoi diagram is created from the 
existing points (ends of the polygon edges) (see Fig.6), it can be observed that this diagram 
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does not meet the geometric constraints, that the polygon edges do not belong to the diagram 

and that the edges are cut. 

In order to introduce the edges of the boundary polygons, it is possible to duplicate 
the existing points. The Voronoi diagram thus has to pass through the edges. However, this is 
5 not yet_sufficient yet-because, although the edges are in the diagram, they still are divided 
(see Fig.7). 

Since the edges of the Voronoi diagram correspond to the mid-perpendiculars of the 
segments defined by two sites, one can try to pas s passing through t hese mid-perpendiculars 
through the ends of the edges defined by the polygons can be tried. For a given end point, the 
1 0 sites of the Voronoi diagram are therefore selected along edges adjacent to this point, 

equidistant therefrom. Since the points are duplicated, one can be certain to obtain an edge of 
the Voronoi diagram on the edge of the polygon. However, if the result obtained is correct in 
the case of a regular polygon (Fig.8a), this is not the case for a non-regular polygon (Fig.8b). 
When the last sites are placed, they are not necessarily at the same distance from the point as 

1 5 the first placed sites. 

Generally speaking, using Voronoi cells for filling the cavity does therefore not seem 

to be suitable. 

SUMMARY OF THE INVENTION 
The method according to the invention allows to simplify simplifies the formation of 
20 a hybrid grid suited to a heterogeneous medium crossed by at least one geometric 

discontinuity of known geometry, in order to form a model representative of fluid flows in 
tnis-a_formation of the medium in accordance with a defined numerical pattern, the structure 
of the medium being known a priori from available data acquired by means of in-situ 
measurements, analyses and/or interpretations of images of the medium (seismic images for 
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example, in the case of a reservoir). 

The method comprises forming at least a first structured grid (irregular CPG type grid 
for example) for gridding of at least part of the medium, forming at least a second structured 
grid for gridding of another part of the medium, forming at least one cavity between the 
5 structured grids including each geometric discontinuity and of sufficient size to allow to form 
allowing forming at least one non-structured grid providing transition between the structured 
grids, and forming each non-structured transition grid by means of power diagrams and by 
imposing conformity of the transition grids with the grid of the walls of each cavity. 

The method is for example applied to a heterogeneous medium where at least one 
1 0 geometric discontinuity is a pipe or a well of known geometry crossing the medium, in 
which case a first structured grid is formed in the medium by respecting considering t he 
discontinuities thereof, a second structured, radial type grid is formed around each well, 
each cavity being delimited defined around each second structured radial grid by 
deactivating grid cells of each first structured grid. 
1 5 The method is also applied for example to a heterogeneous medium where at least one 

geometric discontinuity is a fracture or a fault crossing the medium, in which case a first 
structured grid and a second structured grid are formed in parts of the heterogeneous 
medium, on either side of each fracture, while respecting considering t he discontinuities 
thereof, each cavity formed in order to include a non-structured transition grid therein 
20 being formed by deactivating grid cells of the first and second structured grids, on either 
side of each fault or fracture. 

The use which is made h e r e Usepf the power diagram technique is particularly 
advantageous and convenient in that it allows to appropriately connect appropriate 
connection of t he non-regular structured grids. 
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The invention also relates to a method of simulating, in accordance with a defined 
numerical pattern, the evolution of a process such as fluid flows in a heterogeneous 
medium crossed by at least one geometric discontinuity of known geometry, in accordance 
with a defined numerical pattern, the structure of the medium being known a priori from 
5 available data acquired by means of in-situ measurements, analyses and/or interpretations 
of images of the medium. This simulation is obtained by forming a hybrid grid consisting 
of at least a first structured grid (irregular CPG type grid for example) for gridding of at 
least part of the medium, at least a second structured grid for gridding of another part of 
the medium, at least one cavity between the structured grids including each geometric 
1 0 discontinuity and of sufficient size to allow formation of at least one nea- 

fttructured unstructured grid providing transition between the structured grids, by using 
power diagrams and by imposing conformity of the-transition grids with the-a_grid of the 
walls of said-the cavities, and by solving the numerical pattern in the grid made on the 
medium for modelling simulating t he process. 
1 5 The global hybrid grid is thus obtained by combination of several grid types : a 

structured reservoir grid, a radial grid around each well, also structured, and nen- 
structured unstructured transition grids which connect the previous two grid types. Each 
one of these grids has its own formation methods. 

Using a reservoir simulator of a well-known type, such as ATHOS™ or SCORE™ 
20 for example, for a reservoir provided with a hybrid grid obtained by means of the method, 
allows production simulations to be performed. 



BRIEF DESCRIPTION OF THE DRAWINGS 
Other features and advantages of the method according to the invention will be clear 



from reading the description hereafter of non limitative examples, with reference to the 
accompanying drawings wherein: 

- FigurenFig. 1 shows an example of a structured grid of a faulted reservoir?; 

- Figure- Fig. 2 shows an example of a radial grid around a vertical well?; 

- Figure Fig. 3 shows an example of grid of a horizontal well?; 

- Figure-Fig. 4 shows, in 2.5D, an example of a gridded reservoir where cavities are 
provided for gridded wells, before th e stag e of creation of non-structured interstitial grids 
intended to connect them -the wells t ogether?; 

- Figure-Fig. 5 shows an example of a 2D grid before creation of the transition grid?; 
. pjgttffr-Fig. 6 shows an example of a Voronoi diagram which disregards geometric 
constraints?; 

Pigure-Fig/7 shows an example of a Voronoi diagram whose edges are divided?; 

- Figures -Figs. 8A, 8B show examples of Voronoi diagrams in a regular case and-m 
with any case with conflict respectively?; 

- Figures- Figs. 9 to 14 illustrate the definitions and properties of the power diagrams?; 
r e minded discussed h ereafter, 

Figure Fig. 1 5 shows an example of a power diagram obtained between two 

polygons?; 

- Figure Fig. 16 shows an example of a 2D interstitial grid formed by means of the 
method according the invention, between a structured reservoir grid and structured grids 
formed around a vertical and horizontal well crossing the reservoir?; 

- Figure-Fig. 1 7 shows an example of a global grid obtained with the method?; and 
Figur e Fig. 1 8 shows five wells each p rovided eaeh-with a radial grid, included in a 

gridded 1 5 reservoir, by means nf nnn n tructuredunstructured transition grids based on any 



polyhedral grid cells. 
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DKTATT F,r> DF SI C R 1 PTIQN DET A ILED DESCRIPTION OF THE INVENTION 

Generation of the hybrid reservoir grid is carried out in stages with 
addition/subtraction of elementary grids of different types. 

1) In order to represent the reservoir as a whole, an i, j, k structured grid of a known 
type known to specialists, referred to as CPG, as described in the aforementioned French 
patent FR2,747,490 filed by the applicant assignee , is for example used for e xampl e. 

The reservoir can be faulted with downcreep ofa block in relation to the other. The 
major horizons and faults are first modelled by continuous surfaces from data resulting 
from an interpretation of seismic images of the subsoil or from data obtained during 
drilling (well markers). The geologic structure is then divided into faulted blocks resting 
on these surfaces. These blocks are individually gridded ? and then reassembled. Gridding 
of a block first consists in gridding the edge surfaces, then the inside is populated by 
transftnite interpolation of the edge surface grids. Relaxation techniques are then applied 
to the edge surfaces and to the inside so as to harmonize and to regulate the grid. The grid 
thus obtained strictly respects the considers the h orizons, the faults and it allows to 
rrsprcr.cn t certain representation of unconformities such as pinchouts. It meets all the 
constraints of geologic nature. 

Such a grid is obtained for example in accordance with the method described in the 
aforementioned French p atent 551-2,747,490. 

2) The well trajectories are drawn synthetically or imported. A structured radial 
grid is then generated around each well in order to take account of the particular 
constraints linked with the flows in the vicinity of these-thewells. 
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In the example shown in Fig.2, the structured grid around a vertical well is of circular 
radial type. It is also a CPG type grid. Its generation first consists in sampling a disc at r, 0 
in the horizontal plane. The 2D grid thus obtained is then projected vertically upon the 
various layers of the reservoir grid. Here, the i, j, k of the matrical structure correspond to 
the samplings at r, 9 and z respectively (see Fig. 1 1). 

The grid around a horizontal well (Fig.3) is i, j, k structured, it is of the same type as 
that of the reservoir, except that a well cannot be faulted. It is also obtained by projecting 
vertically upon the various layers of the reservoir grid a 2D grid belonging to a horizontal 
plane. 

3) This radial grid is then inserted around the or around each well in the global 
reservoir grid. A cavity is therefore first created in the reservoir grid by deactivating all the 
grid cells in contact with well grid cells (Figs.4, 5). The space freed between the reservoir 
grid and the well grid must be sufficient to allow convenient formation of a transition grid. 
It can represent for example about the equivalent of two grid cell layers. 

4) A non-structured transition grid is then generated in this cavity (Figs. 16, 1 7) in 
order to connect the structured radial grid around the well to that of the reservoir best 

respecting the constraints linked with the numerical pattern. The user can deactivate the 
grid of a well any time by reactivating the grid cells of the corresponding cavity in the 
reservoir grid. 

Generation of the transition grid uses a technique known in the art, which is for 
example described by: 

- Aurenhammer F., 1987 : Powe r diagrams Diagrams : properti e s Prope rties , algorithms 

Algorithms and applications Applications ; SIAM J. Comput, 16(1), 78-96, 

some notions of which are reminded discussed h ereafter to clarify the description. 
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Weighted point 

Let p be a point of 5R 3 with which a weight w is associated. The pair (p,w) thus 20 
formed is denoted by p and referred to as weighted point. A weighted point can also be 
interpreted as a circle C of centre p and of radius Vw . 

Notion of power 

The power of a point x in relation to weighted point p or to circle C is denoted by 
Power (p , x ) = Power (S , x ) = \px\ 2 — w. \px\ 2 represents here the Euclidean distance between 
p and x. Geometrically, if x is outside the circle and if t is the intersection point of a 
tangent to C passing through x, Power (S)X ) = \xt\ 2 (see Fig.9). 

Power of a point in relation to a circle 

Let C be a circle of eertee- center c and of radius r 5 of equation: 

C(x) = 0whereC(x)xc 2 - r. 

C(x) is the power of point x in relation to C. 

Power diagram 

Let S be a set of weighted points. Each point p of S is associated with a region R(p) 5 
potentially empty, defined by: 

R( p ) = {x e 5R 3 /Power (P)X) < Power^,*), pi G S - {p} } 
This set of regions or cells is called power diagram of S. 
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Radical axis of two circles 

The radical axis of two circles is the line formed by the points having the same power 
in relation to the two circles. In particular, if the two circles have an intersection that is not 
equal to zero, the radical axis is the line passing through the intersection (Fig. 10). 

5 

Use of the power diagrams 

Power diagrams lend themselves well to solution of the problem posed because they 
afford the advantage of having the following properties: the cells that constitute them are 
convex, each cell has a eentre -center (site) and the line connecting the centres-centers of 

10 two adjacent cells is perpendicular to the edge common to the two cells, properties which 
thus satisfy part of the conditions laid down. 

Furthermore, it is advisable that the cells of the power diagram respect the edges of 5 
the boundary polygons. The edges of the polygons therefore have to be edges of the power 
diagram. The constraint imposed for such a diagram is that the edges of the polygons have 

15 to be the edges of a Delaunay triangulation commonly used in the sphere-field of grids. 
This requirement is readily met because the boundary polygons are determined as the 
cavity is created. Once this condition established, it is possible, as explained hereafter, to 
create a power diagram meeting the constraints. 

20 Constrained power diagram : accordance with the given edges 

Let (5) be a set of n edges in the plane belonging to a Delaunay triangulation (Td) (see 

Fig. 11). 

Let (S*) be all the dual edges of (5) in the Voronoi diagram associated with (Td) 15 
(see Fig. 12). 
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Let a s and bj, i = l..n, be the ends of edge S ; . Consider its dual edge S ; * in the Voronoi 
diagram associated with (T d ). It defines a set of circles (Q) of eefitres -centers c ; taken on 
Si* and of radius flc^l = |c, 6, || . Any circle C, of(Ci) therefore comprises &i and b i5 and it 
is an empty circle (Fig. 13). Consider C,' and Cf two circles of (Q). According to all 20 
5 that has been written above, a; and bj are on C J and C ) . a s and b; are thus the intersection 
points of C I and C ) and they are located on their radical axis. In terms of power, it can be 
written in the form: 

C^aO = Cf ( ai ) = C) (bi) = Cf (bi) = 0 

By construction, whatever Q of (Cj), j = 1 ..n, j ± i, the inside of Cj is empty and 
10 therefore contains neither a it nor bi. It is thus possible to write: 
V/,C,.(a,.)>0 
V/,C,(&,.)>0. 

ai therefore belongs to an edge of the power diagram borne by the radical axis of C\ 
and of Cf , bi belongs to an edge of the power diagram borne by the radical axis of C) and of 
15 Cf (Fig. 14). 

Since, by definition, the cells ofa power diagram are convex, wo deduce it is deduced 
therefrom that segment {a b bj] = Si belongs to the power diagram. 

This is true for any Si of (S). A constrained power diagram containing all the 
segments of (S) can thus be constructed. The degree of freedom for constructing such a 
20 diagram is great because the only constraint consists in selecting two eeHfres-craiters_per dual 
edge Si* for any Si of(S). 

It can be noted that, if Si* is reduced to a point, all of (Q) is reduced to a circle of 
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eeatfe- centerc i = Si*. A virtual eentfe- center c iv , infinitely close to the first one, therefore 
to be used in such a case in order to define a second circle. 



Construction of the transition grid by means of the power diagrams 

5 Creation of a power diagram, transition grid, starts with the construction of a 

Delaunay triangulation of the vertex ends of the edges of the boundary polygons. It has 
to be checked that all the initial edges really are edges of the triangulation obtained. The 
Voronoi diagram associated with the vertices of the polygons is then constructed. This 
Voronoi diagram defines, for each edge A of the polygons, a dual edge A'. 
10 Then, for each edge A, two sites of the final power diagram are positioned on A'. 

Since the sites must be weighted points, a weight equivalent to distance \\axf where a is an 
end of A is associated with site x located on A'. The power diagram obtained from these sites 
then respects the edges of the polygons. Figufe-Fig^l6 is an example of a power diagram 
obtained between two polygons with a single site per edge. 
1 5 Once this diagram obtained, all the initial constraints are now met, even if it has to be 

checked that each site is inside its cell. By construction, the sites are positioned on the dual 
edges without specifying in which way. It has so far been shown that a correct power diagram 
is available from the moment that the sites are selected on A', but the possibilities of 
selection of these sites are multiple. It is however possible to guarantee that, if any edge A 
20 intersects its dual edge A', there is at least one power diagram respecting considering the 
boundary polygons and whose sites are inside their cell. 
Results 

Piewe-Fig. 16 shows two examples of 2D grids. The first example corresponds to a 
vertical well situation, whereas the second one shows the case of a horizontal well. It can be 



seen in 
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i these examples that the cells generated for connection of the well and reservoir grids 
meet the constraints stated in the problem : the cells are convex, the sites are inside their cell, 
the edges of the polygons (grid edges) are respected-considered and the line connecting the 
sites of two adjacent cells is perpendicular to the edge common to the two cells. 
5 The 2D stage first allows to generate generation of t he transition grids and thus to 

connect the well grids with the reservoir grid. In order to obtain volume grids, the initial grids 
are converted to 2D grids and a reverse treatment is applied to the 2D results, which allows te 
ebta4n -obtaining a final global 2.5D grid. 

The case where the second structured grids are radial type grids formed around wells 
10 crossing the medium, with delimitation of tefining_cavities around each second grid in order 
to include a transition grid, has been considered. It is however clear that the method is applied 
for gridding of a medium with other types of geometric discontinuities, such as for example 
an underground reservoir crossed by fissures. In such a case, structured grids are formed in 
parts of the medium on either side of the fissures, and cavities are delimited defingdby 
1 5 deactivating grid cells of the two structured grids around each fissure, large enough to form 
non-structured transition grids therein by means of power diagrams and by imposing 
conformity of the transition grids with the grid of the walls of each cavity. 



